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Recent experiments on magnetoresistance in dual spin valves reveal nonlinear features of electronic 
transport. We present a phenomenological description of such nonlinear features (current-dependent 
resistance and giant magnetoresistance) in double spin valves. The model takes into account the 
dependence of bulk/interface resistance and bulk/interface spin asymmetry parameters for the cen- 
tral magnetic layer on charge current, and consequently on spin accumulation. The model accounts 
for recent nonlinear experimental observations. 
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I. INTRODUCTION 

Spin accumulation (spin splitting of the electrochem- 
ical potential) is a nonequilibrium phenomenon which 
is associated with a spatially nonuniform spin asymme- 
try of two spin channels for electronic transport^"— . In 
the simplest case, it appears at the interface between 
ferromagnetic and non-magnetic metals, when current 
has a nonzero component perpendicular to the inter- 
face^. Spin accumulation also appears in more com- 
plex systems, like single or double spin valves exhibiting 
current-perpendicular-to-plane giant magnetoresistance 
(CPP-GMR)^i^ effect, as well as in single or double tun- 
nel junctions. Current-induced spin accumulation is par- 
ticularly pronounced in spin-polarized transport through 
nanoparticles2i^ or quantum dots and molecules^. 

In the case of spin valves based on layered magnetic 
structures, spin accumulation and GMR are usually ac- 
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counted for in terms of the Valet-Fert descriptio 
which the spin accumulation is linear in current, while re- 
sistance and magnetoresistance are independent of cur- 
rent magnitude and current orientation. The descrip- 
tion involves a number of phenomenological parameters 
which usually are taken from CPP-GMR experimental 
data. Originally, it was formulated for coUinear (parallel 
and antiparallel) magnetic configurations, but later was 
extended to describe also current induced spin torque^i 
and CPP-GMR for arbitrary noncollinear geomctrj*i^. 

The Valet-Fert description was successfully applied not 
only to single spin valves, but also to double (dual) spin 
valvesi^, Fl/Nl/Fc/Nr/Fr, where Fc is a magnetically 
free layer separated from two magnetically fixed outer 
layers (Fl and Fr) by nonmagnetic spacers (Nl and 
Nr). An important feature of such structures is an en- 
hanced spin accumulation in the central layer (Fc) for 
antiparallel magnetizations of the outer magnetic lay- 
ers (see Fig.l). Spin accumulation may be then several 
times larger than in the corresponding single spin valves. 
Accordingly, such a magnetic configuration of dual spin 
valves (DSVs) diminishes the critical current needed to 
switch magnetic moment of the central layer, and also 
enhances the current-induced spin dynamic a^'^i^^ . 




FIG. 1: Scheme of a dual spin valve with antiparallel configu- 
ration of outer magnetic layers, Fl and Fr. Fc is the central 
magnetic layer . 



Another interesting consequence of the enhanced spin 
accumulation in the central layer of a dual spin valve 
is the possibility of nonlinear transport effects. Recent 
experimental results^^ indicate that the enhanced spin 
accumulation may cause unusual dependence of magne- 
toresistance on dc current. It has been shown that when 
magnetizations of the outer layers are antiparallel, resis- 
tance of a DSV for one current orientation is lower when 
the Fc layer is magnetized along the Fr one and higher 
when it is aligned along magnetization of the Fl layer, 
while for the opposite current orientation the situation is 
reversed. Moreover, the difference in resistance of both 
coUinear configurations markedly depends on the applied 
current. These observations strongly differ from the pre- 
dictions of the Valet-Fert model^, which gives resistance 
(and magnetoresistance) independent of the current den- 
sity. 

Such a nonlinear behavior may originate from several 
reasons. The Valet-Fert description is based on the as- 
sumption of constant (independent of spin accumula- 
tion and current) basic parameters of the model, like 
bulk/interface resistance, bulk/interface spin asymme- 
try, spin diffusion lengths, etc. This is justified when 
spin accumulation is small and/or change in the density 
of states on the energy scale comparable to spin accu- 
mulation is negligible in the vicinity of the Fermi level. 
Density of states can be then considered constant, i.e. 
independent of energy. Since the density of states deter- 
mines electron scattering rates, one may safely assume 
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that the transport parameters mentioned above are also 
constant. However, when the density of states at the 
Fermi level varies remarkably with energy and spin ac- 
cumulation is sufficiently large, this assumption may not 
be valid, and the parameters mentioned above may de- 
pend on spin accumulatioi*^^. This, in turn, may lead 
to nonlinear effects, like the experimental ones described 
abovoi^. 

The spin accumulation, however, is rather small - of 
the order of 0.1 meV for current density of 10® A/cm^. 
Thus, to account for the experimental observations one 
would need rather large gradient of the density of states 
with respect to energy at the Fermi level. More specif- 
ically, to account the experimental observations, the 
change in density of states should be of the order of 10% 
on the energy scale of 1 meV. Although this is physi- 
cally possible, one cannot exclude other contributions to 
the effect. Spin accumulation can directly change effec- 
tive scattering potential for electrons at the Fermi level. 
Moreover, spin accumulation can also indirectly influence 
transport parameters, for instance via current- induced 
shift of the energy bands due to charging of the layers or 
due to electron correlations, which are neglected in the 
description of the spin accumulation. Since the experi- 
mental results show that the nonlinear effects appear only 
in the antiparallel configuration, where spin accumula- 
tion in the central layer is large, we assume that the indi- 
rect contributions are proportional to spin accumulation 
(at least in the first order). Since, it is not clear which 
contribution is dominant, we present a phenomenologi- 
cal approach, which effectively includes all contributions 
to the observed nonlinear transport. We assume that 
bulk and interfacial resistances as well as spin asymme- 
tries vary with spin accumulation and show that such 
variation leads to effects comparable to experimental ob- 
servationsi^. 

Structure of this paper is as follows. In section II we 
describe the model. Numerical results are presented in 
section III for bulk and interfacial contributions. Section 
IV deals shortly with magnetization dynamics in DSV. 
Finally, we conclude in section V. 



II. MODEL 

Electron scattering rate and its spin asymmetry be- 
come modified when the spin-dependent Fermi levels are 
shifted (eg. due to spin accumulation). All the ef- 
fects leading to this modification can be included in the 
description of charge and spin transport presented in 
Refs.[ll| and[l3, which generalize the Valet- Fert model to 
noncollinear magnetic configurations. We analyze the sit- 
uation when the effect originates from the bulk resistivity 
and bulk spin asymmetry factor /3 of the Fc layer, which 
are assumed to depend on spin accumulation, as well as 
from similar dependence of the corresponding interface 
parameters. Let us begin with the bulk parameters. 

The spin-dependent bulk resistivity of a magnetic layer 



is conveniently written in the form^ 

p^(;) = 2p*(lT/3), (1) 

where p* is determined by the overall bulk resistivity pp 
as p* = pf/(1 — /3^). When the spin accumulation is 
sufficiently large, one should take into account the cor- 
responding variation of p*. In the lowest approximation 
(linear in the spin accumulation) one can write 

P*=Po + <l (.9) , (2) 

where Pq is the corresponding equilibrium (zero-current 
limit) value, and g(x) is spin accumulation in the cen- 
tral layer, which varies with the distance from layer's 
interfaces. To disregard this dependence, we average 
the spin accumulation over the layer thickness (g) = 
(l/d) /p^ g[x)dx. In Eq.([2]) g is a phenomenological pa- 
rameter, which depends on the relevant band structure. 
This parameter effectively includes all effects leading to 
the modification of transport parameters. 
This equation can be rewritten as 

(l + 9*(.9)) : (3) 

where g is a dimcnsionlcss variable related to spin accu- 
mulation, g = (e'^joPo^sf)"^ 5, with jo denoting the par- 
ticle current density and Zgf being the spin-flip length. 
Wc also introduced the dimcnsionlcss current density 
i = I/Iq^ with / = ejo denoting the charge current 
density and /q being a current density scale typical for 
metallic spin valves (/o = lO^A/cm^). The parameter q 
in Eq. ([3]), q = (e/o4f)(7, is a dimensionless phenomeno- 
logical parameter which is independent of current. 

The bulk spin asymmetry parameter (3 becomes mod- 
ified by spin accumulation as well, and this modification 
can be written in a form similar to that in the case of p* , 
i.e. 

/3 = /3o + m , (4) 

where /3o is the corresponding equilibrium value and ^ ef- 
fectively includes all the contributions. When introduc- 
ing the dimensionless spin accumulation defined above, 
one can rewrite Eq. ([¥]) as 

P^pQ + b{g). (5) 

where | = (e/opS^sf )^- 

Similar equations can be written for the interfacial 
resistance R* and interfacial asymmetry parameter 7, 
which define spin-dependent interface resistance as 

i?t(;) =2i?*(lT7)- (6) 

Analogously, we can write the dependence of R* and 7 
on spin accumulation in the form 

R* =Rl+q' g{x,), (7a) 
1 = 10+^' 9{xi), (7b) 
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where g{xi) is spin accumulation at a given interface. 
The constants Rq and 70 are equiUbrium interfacial resis- 
tance and asymmetry parameter, respectively. Relations 
([7]) lead to the following dependence of the interfacial 
parameters on the current density: 

R'' = R*o{l + q't~g{x,)), (8a) 
1 = 10 + ^' i 9{xt) , (8b) 

where q' = (e/oPo^sf) <?' and |' = {eIoPolsi)C- 

The parameters q, ^, q' , and introduced above de- 
scribe deviation from usual behavior of the resistance 
(magnetoresistance) described by the Valet-Fert model. 
These parameters will be considered as independent phe- 
nomenological ones. 



III. NUMERICAL RESULTS 

To calculate resistance and spin accumulation for arbi- 
trary noncollinear magnetic configuration, we apply the 
formalism described in Refs. [ill andfl^ This formalism, 
however, is modified by assuming p* , /3, R* and 7 to de- 
pend on current density (spin accumulation). Therefore, 
for a particular magnetic configuration and for certain 
values of i, q, ^, q', and 5', the spin accumulation has to 
be calculated together with p*, R* , /?, and 7 in a self- 
consistent way. In the first step, we assume equilibrium 
values; p* = Pq and l3 ^ I3q {R* ^ Rq and 7 = 70), and 
calculate the corresponding spin accumulation go{x) in 
the central magnetic layer. Then, we calculate the zero 
approximation of the out-of-equilibrium parameters ac- 
cording to Eqs. (|3]), (O, and/or With these new 
values for p* and /3 (i?* and 7) we calculate the out-of- 
equilibrium spin accumulation in the central layer and 
new out-of-equilibrium values of p* and /3 {R* and 7). 
The iteration process is continued until a stable point is 
reached. Finally, for the obtained values of p*, f3, R*, 7, 
and spin accumulation, we calculate the resistance R of 
the DSV at a given magnetic configuration (see Ref. [H) . 

In all our calculations, magnetizations of the outermost 
layers are assumed to be fixed and antiparallel (like in ex- 
perimenti^). Current is defined as positive for electrons 
flowing from Fr towards Fl. The equilibrium param- 
eters have been taken from the relevant literature (see 
Appendix). 

In this section we apply the above described model 
to two examples of DSV structures. The first one 
is a symmetric DSV with Fl ~ Fr ~ Co(20nm), 
Fc = Py(8nm), and with the magnetic layers sepa- 
rated by lOnm thick Cu spacers. The second consid- 
ered structure is an asymmetric exchange-biased DSV 
similar to that used in experimenlji^, namely Cu- 
Co(6)/Cu(4)/Py(2)/Cu(2)/Co(6)/IrMn(10)-Cu, where 
the numbers in brackets correspond to layer thicknesses 
in nanometers. 



A. Bulk effects 

In this subsection we consider pure bulk effects assum- 
ing q' = and ^' = 0. We start from a symmetric DSV, 
and the corresponding numerical results are shown in 
Fig.[2l First, we analyze the case with q = 0.1 and = 0. 
Figure [DJa) shows how p* varies when magnetization of 
the central layer is rotated in the layer plane. This rota- 
tion is described by the angle 9 between magnetizations 
of the Fl and Fc layers. The higher the current den- 
sity, the more pronounced is the deviation of p* from its 
equilibrium value p^. The current-induced change in Pq 
reaches maxima when magnetic moment of the central 
layer is collinear with those of the outer layers. These 
maxima are different for the two opposite orientations of 
the magnetic moment of Fc layer, as the corresponding 
spin accumulations are different. For 9 = 7r/2, however, 
one finds p* = pg. This is because spin accumulation 
vanishes then due to opposite contributions of both inter- 
faces (for symmetric DSVs). Variation of p* in Fig. [2l^a) 
is shown only for positive current, i > 0. When current 
is negative, the change in p* due to spin accumulation 
changes sign (not shown), as also can be concluded from 
Fig. [Sic). 

The current-induced angular dependence of p* makes 
the resistance of the DSV dependent on the current den- 
sity. As shown in Fig. ^c) , the angular dependence of 
the resistance, becomes asymmetric, i.e. its magnitudes 
in the opposite collinear states {9 = and tt) are dif- 
ferent. Such an asymmetric angular dependence quali- 
tatively differs from that obtained from the Valet-Fert 
description, where the resistance is symmetric. When 
magnetization of the central layer switches (e.g. due to 
an applied magnetic field) from one collinear state to the 
opposite one, one finds a drop (positive or negative) in 
the resistance, defined as AR = R{9 = tt) — R{9 = 0). 
Moreover, when the current direction is reversed, the cor- 
responding drop in resistance also changes sign, as shown 
in Fig [Sic). 

Let us consider now the situation where /3 changes 
with the spin accumulation (current), while p* is con- 
stant, ^ = 0.1 and q = 0. General behavior of f3 and 
of the corresponding resistance with the angle 9 is sim- 
ilar to that discussed above (see Fig[SIb,d)), although 
the sign of the resistance drop for the current of a given 
orientation is now opposite to that obtained in the case 
discussed above, compare Figs[2Ic) and (d). Generally, 
the sign of the drop in resistance may be controlled by 
the parameters ^ and q. 

In real structures, however, both parameters, ^ and 
q, may be different from zero, and the observed be- 
havior results from interplay of the bulk and inter- 
face effects discussed above. To show this, we consider 
now an asymmetric exchange-biased DSV structure, Cu- 
Co(6) /Cu(4) /Py(2) /Cu(2) /Co(6) /IrMn(10)-Cu, similar 
to that studied experimentally. 

Figures [Sja) and (b) show the current-induced angular 
dependence of p* and /? for g = ^ = 0.1. In comparison to 
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FIG. 2: Symmetric dual spin valve Cu- 

Co(20)/Cu(10)/Py(8)/Cu(10)/Co(20)-Cu: (a) angular 
dependence of p* for q = 0.1 and ^ = 0; (b) angular depen- 
dence of /? for ^ = 0.1 and q = 0; (c) angular dependence of 
the resistance (per unit square) for q — 0.1 and ^ = 0; (d) 
angular dependence of the resistance for ^ = 0.1 and q = 0. 
The relative current density i as indicated. 
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FIG. 3: Asymmetric exchange-biased dual spin valve Cu- 
Co(6)/Cu(4)/Py(2)/Cu(2)/Co(6)/IrMn(10)-Cu: angular de- 
pendence of p* (a) and /3 (b) for q — 0.1 and ^ = 0.1, and the 
angular dependence of the corresponding resistance (c); (d) 
dependence of the drops in resistance (per unit square) as a 
function of the reduced current density i for g = 0.1, ^ = 0.1 
[line (1)], q = 0.1, f = 10"^ [line (2)], and q = 10'^, ^ = 0.1 
[line (3)]; (e) drop in the resistance as a function of iq and 
(with reduced current density i); the line covers the points 
where AR = 0. 



the symmetric DSV structure, the difference in the devi- 



ations of both parameters from their equilibrium values 
for = TT and = is now much more pronounced. As 
before, the nonequilibrium values of the parameters cross 
the corresponding equilibrium ones for nearly perpendic- 
ular configuration, 9 ~ 7r/2. The resistance shown in 
Fig.E^c) reveals well defined drop between both collinear 
configurations, and the drop changes sign when the cur- 
rent is reversed. However, the resistance drops arc now 
different in their absolute magnitude due to the asymme- 
try of DSV. 

Figure El^d) shows the resistance drops as a function 
of the current density for three different sets of param- 
eters. For the parameters used in Fig. (3a-c), i.e. for 
g = ^ = 0.1 [line (1)], the absolute value of the drop in- 
creases rather linearly with increasing magnitude of cur- 
rent, although the growth of AR is faster for positive 
than for negative current. In the second case, 5 = 0.1 
and ^ = 10"'^ [line (2)], the dependence remains nearly 
the same, with only a small deviation from the first case. 
Finally, we reduced the parameter q, q ~ 10~^, while 
^ = 0.1 [line (3)]. Now, the dependence strongly differs 
from the first two cases. AR only slightly varies with cur- 
rent and remains rather small. Such a behavior results 
from interplay of the bulk and interface contributions. 
This interplay is presented also in Fig. Elje), where the 
resistance drop is shown as a function of iq and i^. Ad- 
ditionally, the latter figure shows that for any value of q 
there is a certain value of ^ for which A-R = (presented 
by the line). 



B. Interfacial effects 

Now we consider the nonlinear effects due to current- 
dependent interfacial parameters, as given by Eqs.[5] For 
both symmetric and asymmetric spin valves we assume 
that the parameters q' and ^' are equal for both interfaces 
of the central layer. Consider first a symmetric DSV. The 
corresponding results are summarized in Fig. |31 Varia- 
tion of R* , when the central magnetization rotates in the 
layer plane, is shown in Fig. SJa) for g' = 0.1 and ^' = 0. 
The curves below the equilibrium value Rq correspond to 
R* on the left interface, while these above JSq describe 
R* on the right interface. When the central magnetiza- 
tion is close to the collinear orientation {9 = 0,7r), R* 
on the left and right interfaces are significantly different, 
and this difference becomes partly reduced when 9 tends 
to = 7r/2 (for the systems considered). Generally, the 
higher current density, the more pronounced is the shift 
of R* on both interfaces from their equilibrium values. 
The corresponding angular dependence of the DSV re- 
sistance is shown in Fig. for the current densities 
I / Iq = ±3. This angular dependence results in small re- 
sistance drops of opposite signs for opposite currents. As 
shall be shown below, the small value of A-R is due to 
a relatively large thickness of the central layer. Similar 
conclusions can also be drown in the case when g' = 
and only 7 depends on spin accumulation, as shown in 
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Figs.Stb) and (d). 
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FIG. 4: Symmetric dual spin valve Cu- 

Co(20)/Cu(10)/Py(8)/Cu(10)/Co(20)-Cu; (a) angular 
dependence of R* on the left (curves below Rq) and right 
(curves above Rq) interfaces of the central layer for q' =0.1 
and — 0; (b) angular dependence of 7 on the left (curves 
below 70) and right (curves above 70) interfaces of the central 
layer for ^' = 0.1 and q' — 0; (c) angular dependence of 
the resistance (per unit square) for q' = 0.1 and ^' = 0; (d) 
angular dependence of the resistance (per unit square) for 
i' = 0.1 and g = 

For the asymmetric exchange-biased DSVs, we assume 
that both R* and 7 depend on spin accumulation. As 
shown in Fig. [SJa) for g' = ^' = 0.1, there is a relatively 
large drop in resistance for the assumed parameters. This 
resistance drop AR increases rather linearly with the 
current density, as shown in Fig. [5l^b). A small devia- 
tion from the Unear behavior can be observed only for 
larger values of negative current. Calculations for differ- 
ent thicknesses of the central layer, d = 2, 8, and 16 nm, 
show that the slope of the curves representing the resis- 
tance drop as a function of the current density decreases 
as the thickness increases, see Fig.lSJb). In other words, 
the dependence of resistance on current becomes less pro- 
nounced when the central layer is thicker. We note, that 
such behavior was not observed in the case of the bulk 
contribution. This feature arises because for thicker mag- 
netic layers the bulk resistivity dominates the pillar resis- 
tance and suppresses the current-induced effects due to 
interfaces. Additionally, the slope of the curves present- 
ing the resistance drop as a function of current density 
depends on the parameters q' and and can change 
sign for appropriate values of these parameters. This is 
shown in Figs.[5l[c) and (d), where one of the parameters, 
either ^' (c) or q' (d) has been reduced to 10~^. Since q' 
and are of the same sign, their effects are opposite and 
the corresponding contributions may partly compensate 
each other. This is also shown in Fig. ^e), where the 
resistance drop AR is shown as a function of iq' and i^' . 
From this figure also follows that total compensation of 
the contributions to the resistance drop occurs for the 
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FIG. 5: Asymmetric exchange-biased dual spin valve Cu- 
Co(6)/Cu(4)/Py(d)/Cu(2)/Co(6)/IrMn(10)-Cu: (a) angular 
dependence of resistance (per unit square) calculated for cen- 
tral layer thickness d = 2nm, q' — 0.1 and ^' = 0.1; (b) 
dependence of the resistance drop (per unit square) on the 
reduced current density i for q' = 0.1, = 0.1, and for differ- 
ent values of d; (c) resistance drop as a function of the reduced 
current density i for q' — 0.1, ^' — 10^^, and indicated val- 
ues of d; (d) resistance drop vs current density i for ^' — 0.1, 
q — 10"'^, and for indicated values of d; (e) resistance drop 
as a function of iq' and i^' , calculated for d — 2. The line 
covers the points where AR = 0. 



IV. MAGNETIZATION DYNAMICS 

In the analysis presented above magnetization of the 
central layer was in the layer plane. However, when the 
magnetization switches between the two collinear ori- 
entations (due to applied magnetic field), it precesses 
and comes into out-of-planc orientations as well. Such a 
processional motion modifies spin accumulation and also 
DSVs resistance. In this section we describe variation of 
the resistance, when magnetization of the central layer 
is switched by an external magnetic field back and forth. 
To do this we make use of the single-domain approxima- 
tion. We also assume that the magnetic field is applied 
along the easy axis of the central layer, similarly as in 
experiment (sec Fig.[T]). 

Time evolution of the spin moment of central layer is 
described by the Landau-Lifshitz-Gilbert equation 

ds , , ~ . ds 

-T-^-\lg\nosxHeff + asx—. (9) 
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Here s is a unit vector along the spin moment of the 
central layer, 7g is gyromagnetic ratio, /ip vacuum per- 
meability, a is a dimensionless damping parameter, and 
Hcs stands for effective magnetic field. 



dcm 



-ffth, (10) 



which includes external magnetic field (Hcxt) applied 
along e^-axis (see Fig. [1]), anisotropy field (i?ani), and 
demagnetization field (Jfdcm) calculated for a layer of 
thickness d = 2 nm and elliptical shape with the ma- 
jor and minor axes 130 nm and 60 nm, respectively. 
Hth is a stochastic gaussian field with dispersion D = 
(aksT)/ {jgiiQAIsV), which describes thermal fiuctua- 
tions at temperature T, where is the Boltzmann con- 
stant, and V is volume of the central magnetic layer. 
Magnetic moments of the outer layers are assumed to be 
fixed due to much larger coercive fields of these layers. 
Moreover, the torque due to spin-transfer has not been 
included. 
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FIG. 6: Asymmetric exchange-biased dual spin valve Cu- 
Co(6)/Cu(4)/Py(2)/Cu(2)/Co(6)/IrMn(10)-Cu: Minor hys- 
teresis loops of resistance in external magnetic field calcu- 
lated for g' = 0.1, ^' = 0.1, and for diflFerent current densities 
i. Only interface contribution is considered here. 



Figure |6] shows quasistatic minor hysteresis loops of 
the resistance in external magnetic field, calculated for 
asymmetric exchange-biased DSV at T = 70 K. These 
figures are in agreement with the results obtained in the 
preceding section, and also in good agreement with ex- 
perimental observations^. They also show that the drop 
in resistance changes sign when the direction of current is 
reversed. Moreover, one can observe small salient points 
in the hysteresis loops, which appear during the reversal 
process ~ especially at low current densities. These points 
indicate on the minima in resistance at noncollinear con- 
figurations and have been observed experimentally as 
well. 

The minor hysteresis loops appear also in the case 
when the nonlinear effect is due to bulk parameters (not 
shown). Some differences however appear, for instance in 
their dependence on the layer thickness. This suggests, 
that the experimentally observed effects are more likely 
due to interface contribution, which is quite reasonable as 
the spin accumulation is maximal just at the interfaces. 



CONCLUSIONS AND DISCUSSION 



In summary, we have extended the description of spin 
accumulation and magnetotransport in order to account 
for nonlinear magnetotransport in metallic spin valves. 
We assumed the dependence of bulk resistivities, inter- 
face resistances, and bulk/interface asymmetry parame- 
ters on spin accumulation in the central layer. The as- 
sumed phenomenological parameters effectively include 
different contributions leading to modification of the 
spin-dependent density of states at the Fermi levels. The 
obtained numerical results refiect the trends observed ex- 
perimentally. More specifically, the dependence on spin 
accumulation of any of the considered parameters leads 
to an asymmetric modification of spin valve resistance in 
comparison to its equilibrium (zero-current) value. This 
modification results in a drop in resistance when the mag- 
netic moment of the central layer switches between two 
coUincar configurations. Moreover, this drop depends on 
the current density, as has been also shown in exper- 
iment^. Within our phenomenological description we 
can reproduce mainly linear dependence of the current- 
induced resistance drops, with a small deviation from 
the linearity for higher current densities. However, the 
description fails to account strongly nonlinear variation 
of AR, which was observed in some DSV structures at 
high current dcnsitiesi^. To account for this behavior one 
should take into account higher order terms in the ex- 
pansion of the relevant parameters. Additionally, when 
only intcrfacial contribution is taken into account, the 
dependence of AR on current becomes less pronounced 
with increasing thickness of the central layer. Such a be- 
haviour has been observed experimentally^^, too, which 
indicates that the interface contribution to the nonlinear 
effects is more important than the bulk one. 

The resistance drop measured experimentally at the 
current density of / = lO^Acm"^ is about 0.04 ffim^. 
To reach effects of similar magnitude within the inter- 
facial model, as shown in Fig. [5l one needs ^' ^ 1, i.e. 
^ ~ 1.13(meV)~^ (when assuming the effect is due to 
variation of interfacial asymmetry parameter only). If 
direct contribution from spin accumulation would dom- 
inate, then the corresponding change in the density of 
states would be of the order of 10% on the energy scale 
of 1 mcV. This slope may be much smaller in the presence 
of other contributions. 
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Appendix A: Parameters used in numerical 
calculations 

The parameters used in numerical calculations pre- 
sented above have been taken from and ? , and are 
summarized below: 



material p* \[i£l cm] /? i^f [nm] 



Co 


5.1 


0.51 


60 


Py 


16 


0.77 


5.5 


Cu 


0.5 





350 


IrMn 


150 





1 



7 



iiit6rf&CG 


R* 
-ft 


7 






Co/Cu 


0.5 


0.77 


0.542 


0.016 


Py/Cu 


0.5 


0.70 


0.390 


0.012 


Co/IrMn 


0.5 


0.10 






IrMn/Cu 


0.5 


0.70 







TABLE II: Interfacial parameters used in calculations; R* is 
given in fSlm^, and mixing conductance, G^^, in (ffim^)^^. 
Only parameters needed in calculations are given. 
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